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In the standard picture of cosmological structure formation, the Universe we see today 
is evolved under the gravitational instability from tiny random fluctuations. In this talk 
I discuss the onset of non-linearity in the large scale structure formation of the Uni- 
verse when the linear perturbation theory break downs. Using ID Zel'dovich Approxi- 
mation which provides an exact solution for density evolution, I illustrate two effects: 
mode spawning and mode merging and their connection to mode coupling. Those mode 
couplings (quadratic, cubic . . . etc.) from gravitational clustering are in fact what the 
polyspectra (bispectrum, trispectrum . . .etc.) are meant to measure. 

1. Introduction 

One of the main issues in cosmology is quantitative characterization of the 
large-scale structure of the Universe. Recent large sky surveys such as SDSS 
(see http://www.sdss.org and the publications therein) and 2dFGRS (see 
http://www.mso.ctnu.edu.avi/2dFGRS/ and the publications therein) have dis- 
played rich structures of our Universe. In the framework of the inflation paradigm 
1,2,3,4,5^ .j-jjg structure we see today is blown up from tiny inhomogeneities originated 
from quantum fluctuations. This picture is confirmed by cosmic microwave back- 
ground (CMB) measurement of NASA WMAP mission 6 > 7 < 8 . 9 < 10 . The homogeneity 
and isotropy with fluctuations only 1 part in 10 5 of CMB at the epoch of decoupling 
provides the strong evidence that our Universe was very smooth Such initial den- 
sity fluctuations in the simplest inflationary universe constitute a Gaussian random 
field. 



2. Gaussian random fields and Fourier phases 

Gaussian random Fields 11,12 are useful because of its analytical simplicity. One 
particularly interesting property of Gaussian random fields is that the real and 
imaginary part of the Fourier modes are both Gaussian distributed and mutually 
independent, or in other words, the Fourier phases are randomly distributed between 
and 2ir. The statistical properties are then completely specified by its second- 
order statistics: its two-point correlation function £(r) (with zero connected n-point 
correlation functions for n > 2), or alternatively, its power spectrum P(k). 
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In the framework of gravitational instability, a perturbative method can be 
adopted at the early stage of density clustering (see e.g. Peebles 1980, Sahni and 
Coles 1995, Bernardeau et al. 2002). The linear perturbation theory is applicable 
when the density fluctuations are small i.e. its variance of density contrast (<5 2 ) < 1. 
The statistical distribution of the density field such as originally random phases 
remains invariant in the linear regime, except its variance increasing with time. 

The departure of the density field from the linear regime gives rise to phase 
coupling. Second-order statistics, such as power spectrum and two-point correlation 
function, throw away the fine details of the delicate pattern of cosmic structure. 
These details lie in the distribution of Fourier phases to which second-order statis- 
tics are blind. The evident shortcomings of P(k) can be partly ameliorated by 
defining higher-order quantities such as the bispectrum 13 . or correlations of 5(k) 2 
16 . Higher-order correlations and polyspectra find this information term by term so 
that an infinite hierarchy is required for a complete statistical characterization of 
the fluctuation field. 

3. Fourier description 

The mathematical description of an inhomogeneous Universe revolves around the 
dimensionless density contrast, S(x), which is obtained from the spatially-varying 
matter density p(x) via 

= p(x)-po^ 
Pa 

where x is the comoving coordinate, pa is the global mean density 13 . The two-point 
covariance function, which measures the excess probability over Poisson distribution 
between a pair with separation r, is defined as 

£(r) = (5(x)6(x + r)), (2) 

where the mean is taken over all points x. It is also useful to expand the density 
contrast in Fourier series, in which S is treated as a superposition of plane waves: 

5(x) =J2 5 ( k ) c M*k-x). (3) 

The Fourier transform 8(k) is complex and therefore possesses both amplitude \S(k)\ 
and phase (p^ where 

«5(fc) = |<5(fc)|cxp^ fc ). (4) 

The power spectrum is defined as 

(S(k 1 )S(k 2 )) = (2wfP(k)S D (k 1 + fe 2 ), (5) 

which is the Fourier transform of the two-point covariance function via Wicncr- 
Khintchin theorem. We can analogously define the bispectrum as the third-order 
moment in Fourier space: 

(6(k 1 )6(k 2 )S(k 3 )) = (27r) 3 B(k 1 ,k 2 ,k 3 )S D (k 1 +k 2 + fe 3 ), (6) 

which is the Fourier transform of the connected three-point covariance function. 
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4. The linear theory of density perturbations 

Gravitational instabilities is believed to be the driving force of large-scale structure 
formation. When the density perturbation is small, evolution of the density contrast 
can be obtained analytically through the linear perturbation theory from 3 coupled 
partial differential equations. They are the linearized continuity equation, 

85 1 



dt 

the linearized Euler equation 



a Vx-«, (7) 



dv a 1 _ 1 „ , 

-^r + -v = V x p--V x <£, 8 

ot a pa a 

and the linearized Poisson equation 

V x 2 4> = A7rGa 2 p a 5. (9) 

In these equations, a is the expansion factor, p is the pressure, V x denotes 
a derivative with respect to the comoving coordinates x, v = ax is the peculiar 
velocity and <fi(x,t) is the peculiar gravitational potential. From Eq.(7)-(9), and if 
one ignores pressure forces, it is easy to obtain an equation for the evolution of 5: 

S + 2(^)6 -4nGp 5 = 0. (10) 

For a spatially flat universe dominated by pressureless matter, po(t) = l/67rGi 2 and 
Eq.(10) admits two linearly independent power-law solutions 

6(x,t)=b±(t)5 (x), (11) 

where 8o(x) is the initial density distribution, b + (t) oc a(t) oc i 2 / 3 is the growing 
mode and b-(t) oc i -1 is the decaying mode. 

As one can see from Eq.(ll), the growth depends only on time, so the growth 
of density distribution in the linear regime does not alter the phases. 



5. The Zel'dovich Approximation 

In order to examine analytically non-linear effects induced by gravitational cluster- 
ing, we use the Zel'dovich approximation 7 as a clustering scheme, which extrapolates 
the evolution of density perturbations into non-linear regime. This extrapolation 
from the linear theory follows the perturbation in particle trajectories rather than 
in density fields. In the Zel'dovich approximation (ZA), a particle initially placed at 
the Lagrangian coordinate q is perturbed after a time t to an Eulerian coordinate 
x. The displacement of the particle simply depends on the constant velocity the 
particle has when it is kicked off its initial position, and can be written as b(t)u(q), 
so that 



r(q,t) = a(t)x(q,t) = a(t)[q + b(t)u(q)}, 



(12) 
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where r is the resultant physical coordinate, a(t) is the expansion factor and b(t) 
is the growing mode b + (t) from the linear perturbation theory. According to this 
prescription, each particle moves with a constant velocity along a ballistic trajec- 
tory, which resembles Newtonian inertial motion. Note that the peculiar velocity 
according to the ZA is ax = a(t)b(t)u(q) . For an irrotational flow, u(q) can be 
expressed as a gradient of some velocity potential — V$o(g)- 

We focus particularly on the applications of onc-dimcnsional ZA. The ZA in 
ID provides an exact solution of density evolution in that the evolution of planar 
collapse from the ZA has the same solution as from the Poisson equation 18 until 
shell-crossing. The ZA in ID is simplified to 

x ( q ,t)=q + b(t)u(q)=q-b(t)^^-. (13) 

The density contrast can be derived from the conservation of mass pdx = p dq : 



6=*—l=(¥) -1 = 



Po \dq 



!-&(*)- 



dq 2 



1 (14) 
(15) 



dq 

The velocity potential $0(9) can always be disintegrated as 

$0(9) =5^i4iCos(Aj« + Oi). (16) 

i 

6. Mode spawning and mode merging 

Upon using ID ZA, we firstly choose as a toy model the velocity potential a 

$o(<?) = -Ai cos(Aig + an) - A 2 cos(A 2 g + a 2 ). (17) 

When b(t) <C 1, the density contrast can be expressed in terms of only the first 
order from Eq.(15) 

, * - m (^) 

= ai cos(Aig + a\) + a 2 cos(X 2 q + a 2 ) 
^aJ^+aJ^), (18) 



a.\ j \oi 2 

where the round brackets hereafter denote cosine functions 19 in order to show 
clearly the harmonic relationship (see below), a\ — b{£)A\^\ and a 2 = b(t)A 2 X 2 . 
After Fourier transform in Lagrangian coordinate we have only two modes with 
wavenumbers Ai and A2, and phases ot\ and a 2 , respectively. So these are the only 



a We put minus sign for both cosine functions in the velocity potential to make the analytical form 
neat, assuming both A\ and A2 positive, otherwise the phases would have a shift by n. 
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Fourier modes the moment the clustering process begins (which we call the "parent 
modes"). It is clear that the first order does not display phase coupling. When the 
second-order term becomes comparable, 

-*)(#) 

L 2 /2Ai\ 2 /2A 2 \ /Ai+A 2 \ /Ai-A 2 \ 

\2ai J z \2a 2 J \a\+a 2 ) \a\ - a 2 J 

Since the density contrast in ID ZA can be expressed as a power series of d 2 <&o/dq 2 
in Eq. (15), the second-order term reflects the quadratic density fields 20 > 21 . This 
second-order term spawns modes with wavenumbers that are from combinations of 
any 2 wavenumbers from the parent ones, i.e. Fourier modes with wavenumbers 2Ai, 
2A 2 , Ai + A 2 and Ai — A 2 spawned from combination of Ai and A 2 . One important 
feature of quadratic density fields is that the phases of the spawned modes follow 
the same kind of harmonic relationship as the spawned wavenumbers, which we call 
hereafter "wavenumber-phase harmonic relation" . Such relation subsequently forms 
phase associations between Fourier modes and is crucial for bispectrum analysis. 

We can generalize the velocity potential in ID ZA as a sum of cosine functions 
$o = — J2i Ai cos(Ai<7 + on). The density contrast of the first order and that up to 
the second order now become, respectively, 

<5 (1) =£a< cos(A, + a { ) = ]T <h ( ) , (20) 

i 

where ai = b(t)Ai\f, and 

* 1 >+*< a > = E*(£)+E 3 ? 

i jk 

If primordial Gaussianity is assumed, the phases a* of the initial velocity potential 
are uniformly random and independently distributed between and 2ir. Below we 
categorize the effects induced by quadratic density fields. 

• Mode spawning : In Eq.(21) beyond the first order, new modes arc 
spawned from the <5^ 2 ) and correlated phases are created following the 
spawned modes. The wavenumbers of the spawned modes are formed 
from combination of any 2 wavenumbers of the first order (the same 
as the toy model) and the phases follow the wavenumber-phase har- 
monic relation. Some terms may appear even at the earliest stage (when 
b(t)AjAk(XjXk/K) 2 > 2Ai). Those are terms usually involving high- 
frequency modes, i.e. high Xj or Afe. Such phase coupling can manifest 



Xj — Xk 



Xj + X k 



(21) 
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(22) 



itself through phase mapping 22 . For example, a short sequence of modes 
is formed with a constant difference in wavcnumber Afc = A^: 

ajOfc / Xj — Afe A / Aj A ajOfc / Xj + Xk 

where the middle mode is taken from the 5^. Due to the wavenumber- 
phase harmonic relation, this phase sequence has a constant phase difference 
A<j) — otk and can be mapped along a line parallel to the diagonal through 
phase mapping technique. Such coupling between modes with large Afc is 
discussed in Chiang, Coles & Naselsky 22 . 

Mode merging : The newly-spawned modes from the second order are 
not all independent but could merge with modes of the same wavenumbers 
in the first order. Take the following two modes from Eq.(21) as an example, 
one from the 8^ and the other from the S^: 

Aj A ajcik I Xj — At 
,aj' 2 \aj-a ki 



(23) 



If Aj — Afe = Ai, i.e. the wavenumbcr of the newly-spawned mode coincides 
with the one in the first order, these two modes merge to form a new mode: 




where 



-i— - + aia,ja k cos(ctij - afe - aj) ( ) , (24) 

a m = arctan ± ^ ^ = Qfc) / 2 (25) 

di cosaj + djdk cos(aj — afe)/2 

where the amplitude is modulated and the phase is shifted. Such modulation 
in amplitudes and phases proceeds gradually and the consequence is that 
the wavenumber-phase harmonic relation is broken. Thus not all spawned 
modes will enjoy the harmonic relation. Note that in this case the evolution 
can still be at very early stage b{t) <C 1 as long as AjAfe/A, 3> (Ai/AjAk) 1 ^ 2 - 

7. Full treatment of the ID ZA 

To account for the interaction of the parent modes leading to mode merging, we 
can give a full treatment of ID ZA by starting from Eq. (14). We have 

' dx\ 1 



5k = h 



i 



e- tkx dx 



dq / 



(26) 



where u(q) = —d& (q)/dq. We can give the same treatment to Eq. (26) as Eq. (15) 
by expressing it in terms of a power series of b(t): 

* -S^s/** <27 » 
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Following Soda & Suto 23 , we can denote the Fourier transform of u(q) as Uk 

Uk = ^f d<1 U{q)e ~ lkq ' (28) 
so the Fourier modes 5k can be expressed as 

(—ikb) 2 \ 

5 k = -ikbuk H — 2J u k 1 Uk 2 5 D [k - (ki + k 2 )] 

ki ,fc2 

+ t^f J- u kl u k2 Uk 3 5 D [k-(k 1 +k 2 + k 3 )] + --- , (29) 

ki.k 2 ,k 3 

where 5 D denotes Dirac-delta function. For u(q) = —d^o/dq = — AjXj sin(A,,g+ 
atj) from the velocity potential $o = — J2j Aj cos (A jq + otj), 

u k = -iV^V^^A, -k) = -^pLe™* = ^ £ 'W 2+a '> (30) 

j 

and U-k is simply its complex conjugate. 

Equation (29) and (30) provide us the insight into mode spawning from the initial 
potential field $o(<?) and mode merging in gravitational clustering. The Dirac-delta 
functions dictate the spawning of Fourier modes, both the amplitudes and phases. If, 
for any Fourier mode k, there is only one non-zero term in Eq.(29), either from the 
parent modes in the $o(<?) (i-e. from the lst-order term —ikbuk) or any combination 
from parent modes such that there is only one non-zero term in Eq.(29) triggered 
by the non-zero 5 D , its phase ctk shall not change and the amplitude grow steadily 
with b(t) n (according to the order of this combination). If, on the other hand, there 
are more than one non-zero term in Eq.(29), indicating there is more than one set 
of parent-mode combination for the wavenumber k [e.g. Eq. (23)], they will interact 
so that the phase ak will change accordingly. 

Take a toy model as an example, where the initial field has only two modes: 
7 and 11. For the lst-order term of Eq.(29), only and uu exist as the parent 
modes. For the 2nd-order term, those which are non-zero for 5 D are k = 7+7, 7+11, 
11 — 7, 11 + 11; for the 3rd-order term, those which are non-zero are k = 7 + 7 + 7, 
11 + 11 + 11, 7 + 7 - 11, 7 + 7 + 11, 11 + 11 - 7 and 11 + 11 + 7 . . . etc.. The 
corresponding phases of these spawned modes come from multiplication of wj and 
Mn and do not interact with each other. 

However, if we have another toy model where the initial field has two modes, 
7 and 14, again the parent modes are u-j and M14. Now for the Fourier mode 
<5fc=i4, for example, we have not only b(t) M14 from the lst-order term, but also 
b(t) 2 Uj 5 D [k — (7 + 7)] from the 2nd-order term, or even higher terms such as the 
4th-ordcr b(t) 4 v% u 2 14 5 D [k - (14 + 14 - 7 - 7)] and the 5th-order b{tf u% u u 5 D [k - 
(7 + 7+ 7 + 7— 14)]. . . etc.. The evolution for 5k=u is therefore from the vector 
addition of all these terms. Of course the higher order terms come into force and 
dominate at later times. One can then imagine that for an initial density distri- 
bution which has continuous Fourier modes, the modes spawned shall immediately 
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interact with some parent modes, thereby complicating the phase and amplitude 
evolution in Fourier domain. 

8. Conclusion 

In this talk, I have used ID ZA to illustrate some non-linear effects induced by grav- 
itational clustering: mode spawning and mode coupling. The onset of non-linearity 
can be illustrated by expanding the density contrast in terms of power series of b(t). 
The second-order of such expansion reflects quadratic density fields, from which 
modes are spawned with phases following a special wavenumber-phase harmonic re- 
lation. The wavenumber-phase harmonic relation holds as long as the wavenumbcrs 
of parent modes do not interact. In a full treatment of ID ZA, such mode spawning 
and coupling take place not only during the departure from the linear regime, but 
also the entire gravitational clustering process. The complexity of Fourier amplitude 
and phase evolution in gravitational clustering comes from continuous interactions 
between mode merging and coupling with the spawned modes. 
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